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Abstract 



Campisi, Zhan, Talkner and Hanggi have recently proposed [I] the use of the logarithmic oscillator as an ideal 
Hamiltonian thermostat, both in simulations and actual experiments. However, the system exhibits several theo- 
retical drawbacks which must be addressed if this thermostat is to be implemented effectively. 



1 The logarithmic oscillator 

A logarithmic oscillator is a point mass to in a 
central logarithmic potential. The Hamiltonian for 
such a particle is 
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where k B T and b can be considered arbitrary pa- 
rameters for the time being. The Hamiltonian equa- 
tions of motion are therefore 
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This mechanical system has several interesting 
properties. 

In the one-dimensional version of the oscillator, 
it is particularly easy to find the equations of motion 
by direct integration (we will disregard the singular- 
ity in the potential for the moment). From (jl.ip . we 
get the value of the momentum, 
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and using the first of Hamilton's equations of mo- 
tion (O), 
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we get a differential equation which can be solved 
by separation of variables 
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Now, the amplitude of the oscillation is determined 
by the points q a that satisfy the following equation: 
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that is, 
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where j3 represents (k B T) . The period of oscilla- 
tion is just twice the time taken by the particle to 
go from q A to q B . 
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The function in the integral is even, so 
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In the more general case, the motion of the par- 
ticle lies on a plane. If it moves in circular orbits 
around the singularity with a radius r, then its ve- 
locity can be deduced from the fact that the central 
and centrifugal forces must balance, 
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Therefore, the speed 
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does not depend on the radius of the orbit. The 
radius of the orbit is a function of the total energy 
E, because inserting (|1.5[) into (|1 . 1|) , setting q equal 
to r and then solving for r gets us 
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Therefore, the time it takes the particle to complete 
an orbit is 



torb. = = 27T./— — beP E . 
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2 Statistical properties 
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For arbitrary initial conditions, the trajectory 
followed by the oscillator will not usually be a closed 
path, but the particle will never move further out 
than 

for a given energy E, and the time between two 
consecutive maximum distances will be somewhere 
between 2tAB and t or b. (note that both times are of 
the same order of magnitude), 
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2 Statistical properties 

The fact that the speed on a circular orbit does not 
depend on the radius is quite surprising. It implies 
that, if an external perturbation were to relocate 
the oscillator on a new circular orbit, the kinetic 
energy would remain the same and all the energy 
absorbed would be completely converted into po- 
tential energy. 

In a sense, this result can be generalised to the 
oscillator's other trajectories. If we define the virial 
G as 

G = pr, (2.1) 
and calculate its time derivative using (|1.2[) . 
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The time average of the previous formula is 
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and if (dG/dt) t — 0, then the average kinetic energy 
must be 
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whatever the value of El This means that the log- 
arithmic oscillator can absorb an arbitrary amount 
of energy without changing its temperature at all, 
behaving (in a way) like an ideal thermostat. 

Is it true, then, that (dG/dt) t = 0? It certainly 
is, as 
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because G has upper and lower bounds, as one can 
see by noting that G is a continuous function, except 
at the origin. Given that 
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G (r max ) = 0, 

we can infer that G (r) has upper and lower bounds 
in the interval (0, r max ], and (I2.2p is correct. How- 
ever, we must not forget that there is a limiting pro- 
cess involved in (12.31) , and hence it might take a very 
long time for the average kinetic energy to converge 
to ksT/2. In fact, we will argue that this is gener- 
ally the case, and that the logarithmic oscillator is 
therefore a somewhat less-than-ideal thermostat. 

A recent article in the ar^iv [T] argued that weak 
coupling between a system of interest and a loga- 
rithmic oscillator will result in canonical sampling 
of the former's phase space. The dynamics of the 
compound system would then be determined by a 
total Hamiltonian 

H (q, p, r, p r ) = H s (q, p) + H osc . (r, p r ) 
+H int . (q, p, r, p r ) = E, 

where Hg (q, p) is the Hamiltonian for the system 
of interest, H osc (r, p r ) is the one-dimensional ver- 
sion of (jl.ip . and H int is the potential energy of 
the weak interaction between the system and the 
oscillator, which we will assume is negligible com- 
pared to Hs and H osc .. The density of states for 
the logarithmic oscillator is 
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with 5 representing the Dirac delta function. The 
integral turns out to be exactly the same as (|1.4[) . 
so 
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Furthermore, the probability density p for a point 
in the phase space of the system corresponding to 
H s is 

n ~' £ -ff"' , ' )) . (2.5) 

The function (E) represents the density of states 
of the compound system, 

£1 (E) = J 5 (E - H (q, p, r, p r )) dr dp r dq dp. 

(2.6) 

Expressions (|2.4[) and (|2.6I) can be used to convert 
into 
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3 Experiments 
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which is precisely the canonical distribution for Hs ■ 
According to the authors of [T], the logarith- 
mic oscillator thermostat has two obvious advan- 
tages. Firstly, contrary to the popular Nose-Hoover 
thermostat, the dynamical equations of motion are 
Hamiltonian. Secondly, it is possible to design ex- 
perimental setups in which the thermostat is an ac- 
tual physical system. Hoover wrote a reply [2] to 
the first claim arguing that Nose-Hoover mechanics 
are in fact Hamiltonian, and included an example of 
an alternative Hamiltonian thermostat of the Nose- 
Hoover type. Campisi et alii answered explaining 
their claim further in [3j. Here we will be consider- 
ing the second claim instead, that is, we will con- 
centrate on the implementation of the logarithmic 
oscillator as a thermostat, both in experiments and 
simulations. 

3 Experiments 

An experimental thermostat that relies on the dy- 
namics of only a few degrees of freedom is no doubt 
a very interesting system. However, the nature of 
the logarithmic oscillator imposes some serious limi- 
tations which must be taken into account before one 
attempts to design such an experiment. 

The first problem is a consequence of the length- 
scales involved. Assume that we wish to bring a sys- 
tem with N degrees of freedom to the equilibrium 
temperature T. If the kinetic energy per degree of 
freedom is initially off by a fraction a of the energy, 

(£), = U + °>3*« r - 

then the logarithmic oscillator will have to ab- 
sorb at least an amount of energy equal to AE = 
Nak B T '/2. We have seen that the oscillator typi- 
cally covers distances of the order of b exp {(3E osc . }. 
The change in energy implies that the distances cov- 
ered will change by 

Ar max _ — r max {^P l) . (^-1) 

This can be problematic if r max , is initially compa- 
rable to the size of the experimental apparatus and 
the oscillator is cooling the system. 

The enormous changes in lengths imply similar 
changes in time scales. Having assumed a weak in- 
teraction between the system of interest and the os- 
cillator, the effect of the interaction on the latter 
during one period of oscillation should not be sig- 
nificant. The period is 
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where A is a factor that depends on the trajectory, 
but which is of the order of magnitude of \/87r, in 
agreement with fj 1 .6(1 . The change in distances car- 
ries with it a corresponding change in periods of 
oscillation, 

Atp^r. — tper. ll - (^'^) 

Therefore, when the oscillator is cooling down the 
system of interest, it will usually move very far out 
and oscillate very slowly. On the other hand, when 
it is "hotter" than the system, it will squeeze into 
a small neighborhood of the singularity and vibrate 
very quickly. 

Let us illustrate the problem with some num- 
bers. The authors of [1] propose an experiment in 
which a small system composed of neutral atoms is 
contained in a box of length L. The logarithmic os- 
cillator is an ion in a two-dimensional Coulomb field 
generated by a charged wire. 

Assume, for example, that we have a dilute 
gas of 10 atoms of argon at an initial tempera- 
ture To = 3 K and that we wish to bring them to 
T = 1 K. This means that the logarithmic oscillator 
must absorb about 

AE = ^Nk B T Q - ^Nk B T = 30 k B T (3.3) 

units of energy. Let us assume further that the cross 
section of the charged wire has a radius equal to 
10 _3 i. Then the logarithmic oscillator must move 
in orbits with 

' max. 

However, when we insert (I3.3[) into (|3.ip we find that 

Ar max . = r max . (e 30 - l) > 10 10 L. 

If we also take equation (|3.2p into account, it is easy 
to see that we should expect to find the oscillator 
outside the box most of the time. 

4 Simulations 

The wide range of time and length scales affects 
the precision and time of computation of numerical 
simulations as well, but the presence of a singular- 
ity in the logarithmic potential introduces another 
complication in the numerical implementation of the 
oscillator, as stepping over the singularity will usu- 
ally lead to the wrong energy E osc , . 

When the particle is in the vicinity of the sin- 
gularity, the slope dH/dr changes very quickly. If 
the oscillator ends up too close to the singularity, it 
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will feel a great force which will push it away from 
the singularity during the next time step, making it 
skip the area in which the potential would slow it 
down again, unless a very small time step is chosen. 

For the one-dimensional version of the logarith- 
mic oscillator, the problem can be solved by calcu- 
lating the new position of the logarithmic oscillator 
first. If the oscillator has stepped over the singular- 
ity, then expression (|1.3I) can be used to calculate 
the time it would have taken to get to the new po- 
sition, and one can reset its kinetic energy to the 
correct value and calculate the evolution of the sys- 
tem of interest during that time. This solution is far 
from satisfactory, though, because it involves find- 
ing numerical values of the error function every time 
the particle passes the singularity. 

A different approach ((Tj) replaces the logarith- 
mic potential with the approximate potential 

1 fr 2 +b 2 \ 

thereby eliminating the singularity and introducing 
only a slight correction in the density of states for 
low values of E osc .. Unfortunately, this imposes a 
limit on the amount of energy available for exchange 
between the oscillator and the system. If the system 
and oscillator are enclosed in a box of length L, one 
only has about ksT In (L/b) units of energy to play 
with. In order to allow for larger energy ranges, 
one must choose smaller values of b (of the order 
of exp{— 2a3N} if we wish to allow the energy to 
fluctuate by a fraction a either way), and this will 
tend to generate a small neighbourhood of r = in 
which the forces on the oscillator are huge. 

5 Conclusions 

The logarithmic oscillator proposed by Campisi, 
Zhan, Talkner and Hanggi displays very interesting 



properties from the point of view of theoretical sta- 
tistical mechanics. However, before it can be used 
as a thermostat in actual experiments and numeri- 
cal simulations, three problems must be addressed. 
Firstly, the distances covered by the oscillator de- 
pend exponentially on its energy. Given that it must 
not interact strongly with container walls or other 
objects, one would expect that it would be very dif- 
ficult to control such a system in practice. Secondly, 
the vast increase in the period of oscillation when a 
system is being cooled down suggests that the de- 
sired thermostated dynamics will be achieved very 
slowly. Lastly, the presence of a singularity intro- 
duces some technical complications in the numer- 
ical implementation of the dynamical behaviour of 
the oscillator. It seems, therefore, that Nose-Hoover 
dynamics will remain a popular option in molecu- 
lar dynamics at least until the problems mentioned 
here are resolved satisfactorily. 
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